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Abstract 

In recent experiments, a sequence of changes in the wetting state ('wetting 
transitions') has been observed upon increasing the temperature in systems 
consisting of pentane on pure water and of hexane on brine. In this sequential- 
wetting scenario, there occurs a first-order transition from a partial-wetting 
state, in which only a microscopically thin film of adsorbate is present on the 
substrate, to a 'frustrated complete wetting state' characterized by a mesoscop- 
ically, but not yet macroscopically thick wetting film. At higher temperatures, 
one observes a continuous divergence of the film thickness and finally, at the 
critical-wetting temperature, the complete-wetting state, featuring a macro- 
scopic film thickness, is reached. This sequence of two transitions is brought 
about by an interplay of short-range and long-range interactions between sub- 
strate and adsorbate. The critical wetting transition is controlled by the long- 
range forces and is, thus, found by determining where the Hamaker constant, 
as calculated from a Dzyaloshinskii-Lifshitz~Pitaevskii-type theory, changes 
sign. The first-order transition involves both short-range and long-range forces 
and is, therefore, more difficult to locate. While the pentane/water system is 
well understood in this respect by now, a detailed theoretical description of 
the hexane/brine system is hampered by the a priori unknown modification of 
the interactions between substrate and adsorbate upon the addition of salt. In 
this work, we argue that the short-range interaction (contact energy) between 
hexane and pure water remains unchanged due to the formation of a depletion 
layer (a thin 'layer' of pure water which is completely devoid of ions) at the sur- 
face of the electrolyte and that the presence of the salt manifests itself only in 
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a modification of the long-range interaction between substrate and adsorbate. 
In a five-layer calculation considering brine, water, the first layer of adsorbed 
hexane molecules, liquid hexane, and vapor, we determine the new long-range 
interaction of brine with the adsorbate across the water 'layer'. According to 
the recent theory of the excess surface tension of an electrolyte by Levin and 
Flores-Mena, this water 'layer' is of constant, i.e. salt-concentration indepen- 
dent, thickness 5, with 5 being the hydrodynamic radius of the ions in water. 
Once this radius has been determined, the first-order transition temperatures 
can be calculated from the dielectric properties of the five media. Our results 
for these temperatures are in good agreement with the experimental ones. 

PACS: 05.70.Fh; 68.10.-m; 68.10. Cr; 68.45.Gd 
Keywords: Wetting, Surface tension, Electrolytes 



I. INTRODUCTION 

Imagine a volatile substance, the adsorbate, to be at liquid-vapor coexistence and in 
contact with a third phase, the substrate, which might be solid or liquid. There are 
two qualitatively different ways in which the three phases can meet: one possibility, 
known as partial wetting, is that the liquid phase of the adsorbate forms discrete 
droplets on the surface of the substrate, which have non-zero contact angle with this 
surface. In this case, the droplets are connected only by a microscopically thin film 
of adsorbate that covers the substrate. According to Young's equation,*^ the contact 
angle 9 is related to the three interfacial tensions between substrate and liquid (7^/), 
substrate and vapor {■jsv), and liquid and vapor (7^^) and can be calculated from 
Isv — Isi = liv COS 6*. The second possible arrangement of the three phases is that the 
liquid phase forms a macroscopically thick layer on the substrate surface in such a 
way that there is no direct contact of substrate and vapor anymore. This situation 
is known as complete wetting and the three interfacial tensions obey Antonow's rule: 
Isv — Isi = liv From this relation, it follows immediately that the contact angle 6 is 
zero for complete wetting. In his seminal article on this subject, CahrP demonstrated 
the possibility of a transition between the two states, partial and complete wetting, 
as, for example, the temperature is varied. He predicted that, close enough to the 
liquid-vapor critical point of the adsorbate, there will be a wetting transition from 
partial to complete wetting (or drying). His theory also predicts this transition to be 
of first-order nature - by now, it has been shown experimentally that this is usually 
the case for wetting transitions! ^ * ^ * '^ * ^ ^ The discontinuous nature of the transition 
manifests itself in the abrupt jump of the layer thickness at the transition temperature 
T^^i, in the existence of a prewetting line and in the occurrence of hysteresis, one 
of the hallmarks of first-order transitions.'^ 
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Nevertheless, since the early days of the theory of wetting, there had been specu- 
lations about the possibility of a critical (or higher-order) wetting transition, in which 
the wetting-layer thickness diverges continuously to a macroscopic value P^ ^" * ^^ * ^^ * ^ ^*^ 
In this context, it is important to distinguish between long-range and short-range 
critical wetting. For the latter, the film thickness is predicted to diverge logarithmi- 
calljMEl (within mean-field theory), whereas, for the former, a power-law behavior 
according to the relation / ~ AT~^, where / is the film thickness and AT = T^^c — T 
measures the distance from the critical-wetting temperature, is expected P * ^" * ^'^ * ^^ 
While the upper critical dimensionality for long-range critical wetting is less than 
three, which implies that the predictions of a mean-field theory are valid, it is equal 
to three for short-range critical wetting.'^'^^ESl contrast to long-range critical wet- 
ting, the occurrence of short-range critical wetting requires the proximity of a bulk 
critical point .l^^ 

While short-range critical wetting has been seen very recently in methanol/n- 
nonane mixtures,!^ the first experimental observation of critical wetting was for an 
example of the long-range type: Ragil et al. reported a continuous divergence of 
the film thickness for pentane on water.-'^'' From the effective exponent of —1 that 
describes this divergence (see above) and the fact that the location of the transition 
(53° C) coincides with the temperature at which the Hamaker constant changes sign, 
it was concluded that long-range forces bring about the critical transition.'^ In accord 
with this assumption is the considerable distance of the critical wetting temperature 
from any bulk critical point. 

For the system of hexane on pure water, i.e. for a three-layer configuration of the 
form water/liquid hexane/vapor, the critical-wetting temperature was estimated - 
based on calculations of the Hamaker constant - to be T^^c = 96° C and is, therefore, 
too high to be observed using the existing experimental set-up.'^ In order to depress 
the (critical-)wetting temperature, salt (NaCl) was added to the system.'^^E^'^ In- 
terestingly, for salt concentrations of 1.5 mol/L and 2.5 mol/L, two marked changes 
of the film thickness were observed ellipsometrically?^ at low temperatures, the ad- 
sorbed film is only a few A thick; on passing the temperature T^,i, there is an abrupt 
increase of the film thickness to about 100 A. The hysteresis which is observed for this 
transition corroborates that this change of the wetting state is of first-order nature. 
Upon increasing the temperature further, the film thickness grows continuously and 
diverges at Tyj^^^ just as it had been observed for pentane on water. Returning to the 
pentane/water system, it was found that, when heating the system from low tem- 
peratures instead of cooling it down from c, as had been done before, a first-order 
thin-thick transition occurs at 25° C in this system as well.!^ 

On the basis of the conventional Cahn theory, the first-order transition in this 
system was predicted to occur at — 30°Cp^ a modification of the original theory by 
Dobbs,!^ however, which treats the first layer of adsorbate molecules in a lattice-gas 
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approximation and, thereby, fulfills Henry's law, brought the estimate of Ty^^i closer to 
25°C (or, to be more precise, to within the range 13-70°C, depending on the effective 
diameter of a pentane molecule that one assumes. Adopting a value of 4.4 A, which 
follows from the excluded volume used in the Peng-Robinson equation of state, yields 
T^,i = 38°C.).™ 

Within this theoretical setting, the observed sequence of wetting transitions ('se- 
quential wetting') is brought about by the interplay of short-range and long-range 
forces. At low temperatures, both long-range and short-range interactions inhibit the 
formation of a wetting layer. At T^^i (the wetting temperature predicted by Cahn 
theory based on short-range forces alone), the short-range forces start to favor a wet- 
ting layer (in a sense that the gain in free energy from having such a layer outweighs 
the cost of creating an additional liquid-vapor interface), the long-range forces, how- 
ever, still act against the formation of a macroscopically thick film; the result of this 
interplay is a compromise: the mesoscopically thick film, which is present in an in- 
termediate wetting state that has been termed 'frustrated-complete wetting' .I^D At 
T^^c, the effect of long-range forces changes its nature from inhibiting to supporting 
a thick wetting layer and, consequently, the layer thickness diverges. 

This interplay also illustrates why sequential wetting has been observed only for 
such a small number of systems: the window of opportunity for it to occur is relatively 
narrow, and, in most cases, the long-range forces will support wetting before the 
short-range forces do, so that only a first-order wetting transition is observed at the 
temperature at which the short-range forces start to favor wetting (as in standard 
Cahn theory). 

As already mentioned, the location of the critical-wetting transition is relatively 
easily calculated from the variation of the Hamaker constant as a function of tem- 
perature. Within Israelachvili's approximation,!^ only the static dielectric constants, 
the refractive indices, and a typical absorption frequency of the substances involved 
are required. 

The location of the first-order thin-thick transition, however, is more difficult to 
predict: within a Cahn-type theory, knowledge of the so-called contact energy, which 
describes the short-range interaction between substrate and adsorbate (as deduced 
from measurements of the surface pressure, for exampl^^ is needed. Furthermore, 
the conventional Cahn theory underestimates the (first-order) wetting temperature of 
alkanes on water considerabljJ^ - only a modification introduced by Dobb^ allows 
one to predict T^^i (semi-)quantitatively, but it requires an effective diameter of the 
adsorbate molecule on the assumption that the latter is spherical .1^^^ 

For pentane on water, it has been shown that this kind of theory, when combined 
with appropriate expressions for the amplitudes of the long-distance tails of the long- 
range forces,!^ is able to reproduce both wetting-transition temperatures accurately, 
T^,i = 298 K and T^,^ = 326 K.l^ Note, however, that, using the same parameters. 
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the short-range forces alone (Dobbs' theory) would predict T^^i = 311 K, so the long- 
range forces are not quahtatively, but, to some degree, quantitatively important also 
for the first-order transition. 

Another unexpected finding in the experiments on the hexane/brine system, in 
addition to the occurrence of sequential wetting, was that the two transition tem- 
peratures, Tyj^i and T^^c, were shifted in parallel as a function of the salt concen- 
tration.'^ On theoretical grounds, it had been expected that the (critical) wetting 
temperature decreases with increasing salinity of the substrate.'^ This behavior is 
indeed observed in the experiments.'^ The first-order transition temperature, how- 
ever, was expectecP^'^ to remain largely unchanged because the brine/air(vapor) 
and brine/alkane interfacial tensions vary in a similar fashion as a function of the 
salt concentration.!^ It was even hoped that the two lines of transition temperatures 
would meet to form a critical endpoint^'' (like in a simple model system of sequen- 
tial wettin^^ - this expectation, however, was not met: as already mentioned, the 
two lines run in parallel and, at first glance, the two transition temperatures might 
appear coupled.' ^^ * *^" * (In contrast to the situation at T^^c, which is necessarily the 
temperature at which the Hamaker constant W changes sign, there is nothing pe- 
culiar or exceptional about the behavior of W at Ty^^i, a temperature which is not 
determined by the long-range forces alone.'^^'^ In an attempt to illustrate that T^,! 
(or at least ATy^^i, the shift of T^^i as compared to the case of hexane on pure water) 
is determined by the same dielectric properties as T^^c, Bertrand et al. proposed 
a Dzyaloshinskii-Lifshitz-Pitaevskii(DLP)-type theorjJ^ to calculate the change in 
contact energy between substrate and adsorbate on addition of salt for the system of 
hexane on brine.!^ With the aid of a Cahn phase portrait (for hexane on pure water, 
neglecting long-range forces), the contact-energy difference is then converted into a 
shift of the first-order transition temperature. The DLP-type calculation takes into 
account the existence of a depletion layer at the brine/alkane interface, which is al- 
most completely devoid of ions! ^^ * '^'^ -^ The calculation of the free energy per unit area 
therefore involves four layers: brine/ water /liquid alkane/vapor.l^ If the thickness of 
the water layer is assumed to be 2 A, which is consistent with estimates based on the 
Onsager-Samaras theory of the interfacial tension of electrolytes,!^ which attributes 
the ions being repelled from the interface to electrostatic effects due to image charges, 
one obtains first-order wetting temperatures that are in very good agreement with 
the experimental results.!^ Since the thickness 5 does not depend on the Hamaker 
constant, this actually just demonstrates that DLP theory works even on very short 
length scales, for which it was not designed originally, but in no way does it support 
the supposition that the Hamaker constant might determine i as well. 

Within Onsager-Samaras theory, however, 5 is a function of the salt concentration 
CNaCb as Bertrand et al. pointed out.'^^'^ If this dependence is taken into account, the 
agreement of the calculated T^^i with the experimental values deteriorates quickly as 
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CNaCi increases. Bertrand et al. speculatively attribute this behavior to the breakdown 
of the underlying Debye-Hiickel theory for higher concentrationsP^ 

In this paper, we take a somewhat different view and develop a theory that - also 
based on DLP theory - accounts for the modifications of the long-range forces due 
to the presence of salt. The system we consider is similar in spirit to the four- layer 
structure (brine/water /liquid alkane/ vapor); in order to have a complete description 
of our system within the theoretical framework, we also take into account the first 
layer of adsorbed hexane molecules, which is denser than the bulk liquid hexane, as a 
separate, fifth medium and, thus, deal with a five-layer structure (brine/water/dense 
liquid hexane/liquid hexane/vapor). To estimate the thickness of the layer of pure 
water, 5, however, we do not rely on Onsager-Samaras theory, which is not valid in 
the concentration range of interest here (0.5 mol/L < CNaCi ^ 2.5 mol/L), but only up 
to CNaCi ~ 0.15 mol/L. Instead, we adopt an idea that has recently been advanced by 
Levin and Flores-Mena in order to improve upon Onsager-Samaras theory for higher 
salt concentration.'^ They had realized that - in addition to the effect of an electro- 
static repulsion of the ions from the interface - there is a concentration-independent 
depletion layer, the thickness of which corresponds to the radius of the hydrated ion, 
which is also about 5 = 2 A for NaCl.'^^'^ Therefore, the apparent puzzle of having a 
constant thickness of the depletion layer is not due to the breakdown of Debye-Hiickel 
theory (which seems to work fine even for relatively high concentrations of 1 mol/L 
in the theory of Levin and Flores-MeneP^ , but is due to the omission of the effect of 
the hydration- sphere contribution in Onsager-Samaras theory. This theory had been 
designed for the low-concentration regime, in which the depletion-layer thickness is 
dominated by electrostatic effects, whereas it is the hydrodynamic ionic radius that 
sets 5 for the concentration range in which we are interested here. 

Now we argue that the contact energy in the hexane/brine system is still the 
same as in the hexane/pure- water case because, as before, only pure water is in 
direct contact with the adsorbed hexane. What changes upon the addition of salt 
is the long-range interaction between the substrate, now brine, and the adsorbate, 
hexane, across the water 'layer' of constant thickness 5. Based on these insights, we 
will develop a new theory to compute the first-order wetting temperatures. 

The remainder of the paper is organized as follows: in the next section, we will 
explain our approach in detail. The results of our calculations will be presented in 
Sec. III. A discussion of the results and of future tasks in Sec. IV will conclude the 
article. 
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II. METHODOLOGY AND OUTLINE OF THE 



THEORY 

This section contains the theoretical framework by means of which we describe the 
wetting properties of an alkane (hexane in this case) on pure water and on brine, 
respectively. First, we will briefly summarize the equation of state that is used to 
represent and compute the bulk properties of hexane. The fact that hexane and (salt) 
water are hardly miscible allows us to use an equation of state for the pure adsorbate 
and enables us to manage without an equation of state for the substrate. In the 
second subsection, we present the modified Cahn-Landau model for the interfacial 
tensions, which, in turn, determine the wetting properties. In that subsection, we will 
also give a detailed account of how the long-range field is modified by the presence 
of salt and of a depletion layer. For completeness and reproducibility, we also list 
the representative equations for the dielectric properties of all five media involved. A 
third subsection will contain some of the technical details of our calculations. 



A. Equation of state for hexane 

As in several previous works on the wetting properties of alkanes (on aqueous sub- 
strates),!^ we employ the Peng-Robinson equation of state to describe the 
thermophysical bulk properties of hexane.^ According to this equation of state, the 
pressure is given by: 

pRT a{T)p' 

1 - 6p 1 + hp{2 -hpY ^ ' 

where T is the absolute temperature, R the molar gas constant, and p the molar 
density. The excluded volume h is determined from critical parameters (see below), 
while the function a{T) also involves the vapor pressure at the reduced tempera- 
ture t = T/Tc = 0.7. Adopting Peng and Robinson's recommendations, we use the 
following prescriptions for the above-mentioned parameters: 

b = 0.07780 RTjPc, (2) 
a(T) = a{Tc)a(t,uj), (3) 
a(T,) = 0.045724 (/?T,)VPe. (4) 

In the second equation, a; is the acentric factor given by 

a; = -l-logio[P(t = 0.7)/Pj; (5) 
the function a{t,uj) is defined as 

a{t,uj) = [l + K{uj){l-t^/^)f (6) 
K{lj) = 0.37464 + 1.54226CU- 0.26992 cj^. (7) 
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The critical and other required parameters for hexane are* 



Pc 



0.296. 



3.010 MPa 



507.7 K 



(8) 
(9) 

(10) 



B. Model for the interfacial tension 

Here, we distinguish between two cases and present the theories to describe the 
wetting behavior of hexane on pure water and on brine separately. 

1. Hexane on pure water 

For hexane on pure water, we use a model of the interfacial tension that is completely 
analogous to the one employed by Weiss and Widom for pentane on water.l^^ This 
model, in turn, is a combination of Dobbs' modified Cahn theory™ and a treatment of 
the algebraic tails of the long-range forces as proposed by Indekeu et al^ Within this 
model, the interfacial tension is obtained by minimizing the free-energy functional 
with respect to the density profile of the adsorbate, p{z), in conjunction with the 
appropriate boundary conditions (see Sec. II. C for details). The free-energy functional 
reads: 



In this model, the z-axis is perpendicular to the substrate surface and the substrate 
(water) is assumed to occupy the lower half-space, for which z < 0, so that the 
planar water/hexane interface is located at z = 0. For 2; > 0, p{z) denotes the 
spatially varying density of hexane, while p^uik is the density of either bulk phase, 
liquid or vapor. The density of the adsorbate at the substrate surface, p{z = 0), is 
denoted by po- The first term, 70, represents that part of the interfacial tension that 
is due to the self-interaction of the substrate (water /vacuum), but since it is only a 
function of temperature and an additive term contributing equally to each interfacial 
tension (except for 7;„, which does not involve the interaction with a substrate), it 
is irrelevant to the wetting properties of the adsorbate, and no value of 70 needs to 
be specified. In the second term, 0(po) denotes the contact energy, which depends 
only on the surface density of the adsorbate. Within the framework of his modified 
Cahn theory, DobbJ^ determined the contact energy for several ra-alkanes on water 




dz 



00 




(11) 
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from experimental data for the surface pressure, applying a procedure that had been 
proposed by Ragil et al^ using standard Cahn theory. For hexane on pure water, 
we find from Fig. 2 in Dobbs' paper: 

0(po) = [-4.16P0& + 2.15(po&)'] PA. (12) 

Here, A is a characteristic length scale of the system (approximately 34 A) defined 
by A = {c/{b'^Pc)Y^'^, where c is the influence parameter (or simply the coefficient of 
the square-gradient term). By matching to the experimental surface tension data of 
n-alkanes, Carey et alW found that c can be estimated from the parameters a and 
b in the Peng-Robinson equation of state for n-alkanes of short and medium chain 
length (on the basis of a corresponding-states idea). Their relation reads^ 

c = 0.27A^^^/^a62/3 + 7.25 x IQ-^" J m^ mol"^ (13) 

with A^^ being Avogadro's constant. 

The third term in Eq. is the continuum part from standard Cahn-Landau 
theory, where A/ measures the excess free energy corresponding to the local density 
over that of either bulk phase. It is given by A/(p, pbuik) = f{p) - PIJ'buik + Pbuik] here, 
f = F/V is the Helmholtz free-energy density, p^tuik the bulk chemical potential, and p 
the local density of hexane, p{z). The coefficient c is the influence parameter as given 
by Eq. (|T3|) . The integration is to be performed from Az to infinity in this case, where 
Az is the 'thickness' of the first layer of adsorbed hexane molecules.'^ Treating the 
hexane molecule as spherical and estimating its diameter a from the excluded-volume 
term in the Peng-Robinson equation of state, one obtains Az = a = 4.4 The 
first layer of hexane molecules {0 < z < Az) is considered explicitly in a lattice-gas 
approach and contributes the fourth term to Eq. (|TT|l . which is a discretized version 
of the third. The density at a distance Az from the substrate is denoted by pi. 

The last term contains the long-range field, which - in a first approximation - 
couples linearly to the density. Following Indekeu et al., we adopt a lower cutoff 
of Zc = 0.2A, which worked well for pentanep^'^ for hexane as the adsorbate, we, 
therefore, have Zc ~ 6.8 A. This value is also very close to 1.5 a ~ 6.6 A, the distance 
of closest approach of adsorbate particles that are not in the first layer, but part of 
the continuum in the modified Cahn theory. The amplitudes of the first couple of 
leading terms in a l/^r-expansion of the long-range forces, 03 and 04, are calculated 
from DLP theory as follows: invoking Israelachvili's approximation,!^ which amounts 
to the assumption that all media involved have the same characteristic absorption 
frequency Ue and that the dielectric spectrum of substance j can be represented by 
ej[iQ = 1 + (n| — 1)/(1 + (^/cue)^), with rij being the refractive index and C denoting 
the frequency, the Hamaker constant of a water(l)/liquid hexane(3)/vapor(2) system 
is given byP^ 

3 / 53(0)- 5,(0) \ / 53(0)- 5^(0) \ 

4 ^ \e.,{Q)+eM) 1^3(0) + 62(0) j 
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3hu!e (n| — nl) {n 



nn 



8V2 (nl + nl)'^' {nl + nj)'^^ Unl + nf + {nl + n^)'^' 



(14) 



In this equation, ks is Boltzmann's constant, while h denotes Planck's constant. The 
coefficient 03 is related to W hj = —W/{6ttpi), where the denominator should 
actually contain the difference pi — p^, but in view of the large distance from the bulk 
critical point of hexane, py can safely be neglected.!^ 

To calculate 04, we follow the approach of Bertrand et alW^ and consider a 
four-layer structure water (1) /dense hquid hexane(4)/liquid hexane(3)/vapor(2) be- 
cause, right at the water/hexane interface, there is a thin layer (of approximately 
one molecular diameter thickness) of hexane whose density is higher than the bulk 
liquid density at the respective temperature; cf. Fig. 1 (a). Since the layer of liquid 
hexane, the thickness of which we denote by /, is assumed to be much thicker than 
just one molecular diameter a in the frustrated-complete wetting state (and certainly 
so in the complete- wetting state), we expand the full free energy per unit area of 
the four-layer structure, as given by Mahanty and NinhanJ^ and by Parsegian and 
Ninham,^ in powers of a /I and truncate the resulting series after the linear term, 
since a ^ I. The Hamaker constant and, therefore, the expression for 03 remain 
unchanged by this consideration: the l~'^-teim. gives the Hamaker constant just as in 
Eq. fll4j) . while the cx/Z^-term results in a coefficient B as appearing in the following 
expansion of the long-range part of the free energy per unit area: 

n ^ W Ba , , 

^^^^^'") = -12;^+12;^' ^'^^ 
where our result for B is given by: 

^ _ ^ /^3(o)h(o)-£4(o)]h(o)+54(o)]^^ ^^3(0) -^^(o)' 




£4(0)[£l(0)+£3(0)]' ; V^3(0) +52(0)7 

■^^{nl-nl){ni-nl)A, (16) 



with A being: 
A 





{nl- 
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[nl - nl) 


2{nl 


- nl) {nl 




3/2 (^2 + ^2 _ 2nl) 
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) {nl - nl) 
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rti ( 


nj + nl- 


2nl) 


^ {nl + nl — 2nl) 
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{nl - nl 


)inl 


-nl-nl + nl) 



{nl - nl)"" {nl + nj)''^ {nl + n§ - 2rz|) 
— + nj {nl + nl + nf) — n\ (ri| + 2nlnl) 

{nf — nl)'^ {nj + nl)^^'^ {nj + nl — 2nl)'^ 
-nl {nj - Snlnl + 2n|) + nl {nj - 2nlnl + 2^4) 
{nf — nl)^ {nl + JT-I)^^^ {n\ + ^3 — 277-4)^ 



(17) 
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We remark that, since B is approximately of the same order of magnitude as W 
(10"^'' J) far from T^^c (it is actually four times larger than W at T = 0°C and ten 
times larger than W at T^, i), the second term in Eq. (jlSp is indeed typically smaller 
than the first by a factor oia/l. In contrast, near T^^c, W is negligible and B becomes 
all- important. 

The amplitude 04 is related to B via 04 = Ba /{Aixpi), where again, has been 
neglected in the denominator since, far from the critical point of hexane, pi ^ Pv Ta- 
ble I contains the representative equations for the dielectric properties involved in the 
above expressions. The static dielectric constant and the refractive index of the hex- 
ane layer near the water surface is estimated using the Clausius-Mossotti and Lorenz- 
Lorentz equations, respectively, assuming that the density of hexane is enhanced by 
about 12% compared to the bulk density of the liquid.!^^^ Numerically, the results for 
B from the above expression are very close to those obtained from the corresponding 
relation derived by Bertrand et al.^^ which, in addition to the truncated a/l- 
expansion, assumes £3 ~ £4 as well as ~ n4 and is, therefore, limited to quadratic 
order in the terms Ajk{i(), where Ajk{iC) = {£j{iC) ~ ^k{K)) I (^j(^C) + ^fe(^C))- This 
additional approximation does not allow one to recover the formally correct result 
for the Hamaker constant in this four-layer calculation. Bertrand's approximate re- 
latiorJ^ reads: 

5 ^ ^ A;bTA32(0)A4i(0) + ^ ^a;eiV(3, 2, 4, 1), (18) 



where 

W,^?,J,A;) 







H- 


-<) 








1/2 


[n] + 


1/2 




{n} + nl) 


+ [n] + nl) 


1/2" 



(19) 



In the following, we propose an alternative approach, which will prove advan- 
tageous in the next subsection concerning the case of hexane on brine. It, too, is a 
quadratic approximation in the quantities Ajk but recovers the exact limits for cr — 
and a 00. According to Mahanty and Ninham,'^the expression for A3.4i((t) in a 
four-layer structure, in which medium 3 is of fixed thickness and media 1 and 3 are 
separated by an intervening layer of medium 4, whose thickness is a (cf. Fig. 1 (a)), 
reads (omitting the frequency-dependence of A for brevity) 

, , A34 + A4ie-''^/' , , 

A3.41 = 7 . 77, 20 

J 1 + A34A4ie— ^/'' ^ ^ 

where x is a dimensionless variable of integration. Note that A3.4i(0) = A31 and 
A3 ;4i(oo) = A34. In Bertrand's approximation,!^ which is conform with the expansion 
proposed by Mahanty and Ninham,^ this expression reduces to A3;4i(cr) ^ A34 + 
A4ie~°'^/', which yields the exact result for a 00, but not for a ^ 0. In the limit 
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0" 0, one has A3.4i(0) ~ A34 + A41, which is not exact but correct to order A^. We 
now observe that A41 = A31 — A34 + O(A^), using Eq. at a = 0. Replacing A41 
by this expression in Eq. (plj) . we obtain A3.4i(cr) ~ A34 + (A31 — A34) e~'^^/\ which 
is correct to order A^ and yields the exact results in both limits a — > and a — > 00. 
It, therefore, allows one to obtain the exact Hamaker constant W in addition to an 
expression for B within an approximation to quadratic order in Aj^, in the context 
of an expansion to first order in a /I. The latter expression reads: 

5 ^ ^ A;BrA32(0) [A3i(0) - A34(0)] + ^ hu, [iV(3, 2, 3, 1) - iV(3, 2, 3, 4)] . (21) 

Numerically, Eqs. ()18|1 and 1)211) give very similar results; to within a resolution of 0.5 
K, there is no detectable difference in the resulting T^,i. 

In the actual process of minimizing the free-energy functional Eq. (|TT|) . we use 
representative equations for the amplitudes 03 and 04, which were obtained by linear 
regressions to the results of Eqs. (fT^ and (j2II), and which, for hexane on pure water, 
read: 

03 = (4.68I X 10-2^J K-^(T - 273.15 K) - 4.460 x 10~2^j) / (67rp,) (22) 

04 = (-5.432 X 10-32 Jm K-^{T - 273.15 K) + 8.987 x lO^^^Jm) / {A-npi) . (23) 

Using the above expressions, we obtain T^, 1 = 346 K and T^^c = 369 K in nearly 
perfect agreement with the values that were extrapolated from the experimental tran- 
sition temperatures for hexane on brine using different concentrations of salt (2.5, 1.5, 
and 0.5 mol/L) down to CNaCi = mol/L.'^ 

2. Hexane on brine 

As already mentioned in the introduction, the fact that ions dissolved in water carry 
a hydration sphere around them, which prevents the ionic centers from approaching 
the (salt) water/alkane interface any closer than the distance set by the radius of 
this sphere, creates a depletion 'layer' of pure water near this interface.^ Despite the 
fact that DLP theory was derived for much larger separations (thicker films) than the 
molecular dimensions we are dealing with here, it seems to work reliably even in these 
extreme cases of very short length scales. With this additional 'layer' of pure water 
present, we now consider a five-layer structure consisting of brine(l) /water (4) /dense 
liquid hexane(5)/liquid hexane(3)/vapor(2); see Fig. 1 (b). In this calculation the 
brine and vapor phases are semi-infinite slabs, while the water layer is of thickness 
5. The layer of dense liquid hexane, again, has molecular dimensions set equal to the 
diameter a of a hexane molecule, and the film of liquid hexane is of thickness /. 

The origin of the z-axis, which marks the border between substrate and adsorbate, 
is now located between the layer of pure water and the molecular layer of dense liquid 
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hexane. Therefore, the major formal difference between layers (4) and (5) is that the 
water layer is part of the substrate {z < 0, but this part is not considered explicitly in 
the free-energy functional), while the first layer of adsorbed hexane molecules clearly 
is part of the adsorbate (2; > 0). In the derivation of Eq. (fTHIl . it was assumed that / 
is much larger than a, so, strictly speaking, this approximation is applicable to the 
frustrated-complete and complete wetting states, but not to partial wetting. In this 
latter state, however, the cutoff Zc being larger than a ensures that the substrate- 
adsorbate (dense liquid hexane, z < a) interaction is entirely accounted for by the 
contact energy. 

In the present five-layer calculation, the long-range interaction of the substrate 
(brine) with the adsorbate, not the one between the water layer and the adsorbate, 
however, now involves the additional distance 5, the thickness of the intervening water 
layer. The complete expression for A3.451 [S, a) with two intervening layers 4 and 5 of 
thickness 6 and a, respectively, now readsP^ 



A35 + A5;4l(^)e— 
1 + A35A5.4i(())e '^^Z' 



where A5;4i((5) is given by 



A54 + A4ie-^-/' 
^'■'''^^^ = 1 + A,4A4ie-^^Vr (25) 



Analogously to what was done for A3.4i(cr) in the previous section, A^.4i{6) is now 
approximated by A5.4i(5) ~ A54 + (A51 — A54)e~^^/^, which is correct to order A^ 
and ensures proper behavior in the limits 5 — > and S ^ 00. Thus, we have 



A35 + 



A54 + (A51 - A54) e 



-5x/l 



^ — crx/l 



A3;45i(5, cr) ^ , , ^ — T A X .-^wn (26) 



1 + A35 [A54 + (A51 - A54) e-^-/'] e---/' ■ 

Taking the limits 6-^0 and cr ^ at the same time leads us to the identity 
A3;45i(0, 0) = (A35 + A51) / (1 + A35A51) = A31, from which we obtain - to quadratic 
order in Ajk - the approximation A51 ^ A31 — A35. Therefore, A3;45i(5, a) now, again 
to quadratic order, becomes 

A3;45i(5, a) ^ A35 + A54 e-'^-'/' + (A31 - A35 - A54) e-('^+'^)-'/'. (27) 

Note that this is exact in the two limits cr — 00 and a = S = 0. From here, 
we proceed by considering the limit 6 — > 00, for which A3;45i(5, a) must transform 
into A3;54(cr). The resulting approximation for the remaining four- layer structure 
(4/5/3/2) is A3.54(ct) ^ A35 + A54e-'^^/^ which, in the limit a ^ 0, yields A54 ^ 
A34 — A35. Substituting this into Eq. (j27j), we arrive at 

A3;45i(5, cr) ^ A35 + (A34 - A35) e-'^^/' + (A31 - A34) e-^'^+^^^Z^ (28) 
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which is the only approximation correct to quadratic order in Aj^ that describes all 
limits 5 — ^ 0, oo and cr — 0, cxd exactly. Incidentally, note that this approximation 
leads to expressions for the free energy per unit area in which the contributions 
proportional to A32A35 /"^ A32 (A34 - A35) {l+a)-^ and A32 (A31 - A34) {l+a+6)-^ 
are easier to interpret than in the textbook expressions, such as, e.g. Eq. (5.8) in the 
monograph by Mahanty and Ninham.l^ Within this approximation scheme, the free 
energy per unit area arising from long-range interactions is 



k T °° 
7LR(/;a,5) = -^E' 

n=0 



A32A35 + -2 A32 (A34 - A35) 

[l + a/l) 



+77^4x777^^32 (A31 - A34) 

{l + a/l + d/l) 



(29) 



where the Aj^ depend on imaginary frequencies i(n, with (n given by Cn = 27inkBT/h, 
and the contributions of these frequencies are to be summed up.^ The prime on the 
summation symbol indicates that the term corresponding to n = is to be multiplied 
by a factor of 1/2. Thus, the presence of two intervening layers gives rise to a free 
energy per unit area that involves contributions proportional to l/(/ + 5)^ and to 
+ 6 + cr)^, respectively. Now, unlike for a /I, we will not expand into powers 
of 6/1 and truncate after the linear term (since this approximation would be valid 
only for frustrated-complete wetting and for complete wetting) in order to keep the 
following physically plausible picture: while the terms involving 03 and 04 concern 
the interaction of the adsorbate with the topmost layer of the substrate (pure water) , 
the tails of the long-range field exerted on the adsorbate by brine will be of the form 
ay{z + (5)3 and ay{z + 6^. Expanding the above equation to linear order in terms 
of a// yields 



n=0 



A32A34 ]~^^'^ (A34 — A35) 



^ rA32 (A31 - A34) - 77-^^— 3 A32 (A31 - A34) 



(30) 



[i + 6/iy i{i + 6/iy 

W Ba W B'a 

(31) 



127r/2 127rP 127r(/ + 5)2 l2T^{l + 5f' 

This equation is seen to consist of four contributions: the first and the second term 
result in the familiar expressions for W and 5, respectively, as given by Eqs. (jl4j) 
and (PT|) . The third and fourth terms, involving W and B\ however, act over a 
distance (/ + 5) and, thus, represent a qualitatively new contribution to the long- 
range interaction free energy. In sum, we obtain the following free-energy functional, 
which, after minimization, will give the interfacial tensions: 



7[P] = 7o + 0(po) + J|A/(p,pb„ife) + ^ (^^^ jdz 
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|a/(po, pbuik) + I [(po - pi)/Azf^ Az 



oo 



- + - j piz)dz - j (^^-^ + j p(.)rf. 

-^PoAz. (32) 

Analogously to the definition of the unprimed coefficients, the primed ones are given 
by 03 = —W'/Qupi and a'^ = E'cr/Anpi, so we now have: 

W ^ ^A;BTA32(0)A34(0) + ^^a;eiV(3,2,3,4), (33) 

B ^ ^A;bTA32(0)[A34(0)- A35(0)] 

+ ^ huj, [iV(3, 2, 3, 4) - iV(3, 2, 3, 5)] , (34) 

W ^ ^A;bTA32(0)[A3i(0)- A34(0)] 

+ ^ huj, [iV(3, 2, 3, 1) - iV(3, 2, 3, 4)] , (35) 
B' = 2W'. (36) 

Note that for the critical transition, from the frustrated-complete wetting state to 
complete wetting, / is very much larger than S and the Hamaker constant, which 
changes sign at this transition, is given by + W. In Table II we have compiled 
representative equations for Og and function of temperature for a few selected 

concentrations of salt. For the primed terms, involving a'^ and a'^, respectively, we 
use the same cutoff Zc as for the unprimed terms. The actual value of Zc is very close 
to 1.5(7, the distance of closest approach of particles in the region which is treated 
as a continuum in the modified Cahn theory employed here. To account for the 
interaction between brine and the first layer of adsorbed hexane molecules across the 
water layer of thickness S, we include the last term in Eq. ()32j) . Therefore, the only 
unknown parameter introduced by extending the theory from having brine instead of 
water as a substrate is 6. It is this last term in Eq. (jH^ that, numerically, accounts 
for the major difference between having water and having brine as the substrate. The 
sensitivity of our theory to the actual value of 6, which will be discussed in detail in 
the next section, can be traced back to this term. 



C. Technical details 

Given the appropriate free-energy (per unit area) functional, Eq. (fTTj) or Eq. (jH^ . 
respectively, we proceed to calculate the different interfacial tensions relevant to the 
wetting behavior. For the critical transition, the three interfacial tensions between 
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substrate (water or brine) and liquid hexane (7s/), substrate and hexane vapor {■jsv), 
and between liquid hexane and the vapor phase (7^) matter. According to Young's 
equation, complete wetting is reached once 7^^, = 7^; + 7;^ (Antonow's rule). In a 
sequential wetting scenario, however, the surface free energy per unit area for the 
frustrated-complete wetting state, 7^^, will be extremely close to the sum 7^; + jiv^ 
so that it becomes very difficult to locate the critical wetting transition exactly via 
this route. Fortunately, for long-range critical wetting, Tw,c is simply determined by 
the temperature at which the Hamaker constant changes sign, which allows one to 
calculate T^^c very easily without any numerical minimization. For the first-order 
transition, we have to compute the free energies of an sf-interface with a mesoscop- 
ically thick film and of an sf-interface without such a layer. Only one of these two 
configurations will be stable, the other one just metastable, except right at T^,i, 
where both are equal in free energy per unit area. 

To obtain the four different interfacial tensions, we prescribe simple initial trial 
profiles that resemble the final structure and then minimize the free energy of the sys- 
tem under the respective boundary conditions (see below) using a conjugate-gradient 
method. All calculations were done using 2000 points distributed evenly over a dis- 
tance of 200 A. Several initial thicknesses of the mesoscopically thick film (60-150 A) 
were tried, but no difference with respect to T^^i could be detected. 

The boundary conditions are p(0) = pi,p{oo) = for the 'free' hquid-vapor 
interface and p{oo) = pi or p(cxo) = p^ for the si and su-cases, respectively, while p(0) 
is to be obtained in the process of minimizing the free energy of the whole system in 
these cases (natural boundary conditiorP^. 

Note that the long-range terms only contribute significantly to 7s„ if / > Zc-, i.e. for 
the frustrated-complete and complete wetting states, but not for partial wetting. The 
long-range interaction between brine and the first-layer of adsorbed alkane molecules 
(last term in Eq. ()32|) ) is very important and, therefore, taken into account in all 
wetting states. 

III. RESULTS 

The main results of our calculations are shown in Fig. 2, which contains experi- 
mentaP^ and theoretical first-order and critical-wetting transition temperatures for 
hexane on brine as a function of the salt concentration. The open circles are the two 
critical-transition temperatures as determined experimentally for CNaCi = 1-5 mol/L 
and CNaCi = 2.5 mol/L. The solid line represents the loci where the Hamaker con- 
stant, as computed from Eq. (fT^ using the equations listed in Table I, changes sign. 
The agreement with the experimental data is very good.*^ The open diamonds de- 
note the three first-order transition temperatures that were obtained experimentally 
for CNaCi = 0.5, 1.5, and 2.5 mol/L.l^ Linear extrapolation down to CNaCi = mol/L 
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yields a first-order transition temperature of 73° C for liexane on pure water The 
filled squares (connected by a dashed line as a guide to the eye) are the theoretical 
values of T^^i as computed from the theory outlined in Sec. II. B. 2 using a thickness 
of the depletion layer of 6 = 1.9 A. As can be seen in Fig. 2, the agreement with 
experiments is as good as for the critical-wetting transition temperatures, but a slight 
bend in the theoretical curve, which leads to somewhat larger deviations at higher 
salt concentrations, is visible. For CNaCi = mol/L, the first-order transition temper- 
ature is calculated according to the specifications given in Sec. II.B.l. In this case, 
the agreement with the extrapolated 'experimental' value of T^,! = 346 K is perfect. 

Naturally, the question of the thickness of the depletion layer does not arise for 
CNaCi = mol/L; for CNaCi > mol/L, however, it does. In their theory of the 
increase of the surface tension of water on addition of a strong electrolyte, like NaCl, 
for example. Levin and Flores-Mend^ used the value of 6 = 2.125 A, which led 
to a good description of the experimental surface-tension data given by Matubayasi 
et a/.E3 for aqueous solutions of this salt. The theory of Levin and Flores-Mena, 
however, is quite insensitive to the actual value of 6: as long as 6 is of the order of 2 
A, it describes the experimental data quite well, as can be seen in Fig. 3 (in which 5 
ranges from 1.75 to 2.5 A). 

In sharp contrast to this behavior, the first-order transition temperature is rather 
sensitive to small changes of 6. In Fig. 4, we have compiled the values of T^,i(cNaCi) for 
6 = 1.75, 1.9, 2.125, and 2.5 A. While these thicknesses are still suitable for describing 
the surface tension of NaCl solutions (cf. Fig. 3), the predictions of T^^i are not 
satisfactory, except for 6 ~ 1.9 A. (Note, however, that because of the curvature of 
^u;,i(cNaCi), choosiug 6 = 1.75 A yields good agreement at higher salt concentrations, 
e.g. for CNaCi = 2.5 mol/L.) In fact, if one chooses 6 large enough, one may even create 
- purely hypothetically, of course - the elusive critical endpoint, where T^^i coincides 
with T^^c- 

For comparison, we also compute 6 from Onsager-Samaras theory as Bertrand et 
al. didpSEa the prediction of i(cNaCi) within our approach (see open diamonds in 
Fig. 4) shows the same downwards-bend as in their theory. The disagreement of cal- 
culated and experimental values of i is therefore mainly due to the inapplicability 
of Onsager-Samaras theory, which predicts a varying 5 because it neglects the hydra- 
tion sphere and bases the estimate of 5 solely on electrostatic considerations. Doing 
so is correct for very low concentrations of salt, but certainly not for the concentra- 
tions present in the experiments on sequential wetting. In conclusion, what causes 
the disagreement is not the breakdown of Debye-Hiickel theory, as Bertrand et al. 
proposed,!^ even though this theory is definitely not applicable for the hulk proper- 
ties of a salt solution in this concentration range. As Levin and Flores-Mena aigue^ 
the reason why a theory of the surface tension based on Debye-Hiickel theory still 
works for these relatively high salt concentrations is that, in the canonical-ensemble 
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approach introduced by LevinP^ the increase of the interfacial tension is given by 
the difference between the bulk free energy of a homogeneous system and that of 
an inhomogeneous system which features the hquid-vapor(air) interface. Therefore, 
most bulk effects simply cancel. 

IV. DISCUSSION AND CONCLUSION 

In the preceding sections, we have demonstrated that a DLP-style theorj^^ combined 
with Israelachvili's approximation^^ and a modern account of the surface tension of 
electrolyte^^ is able to describe not only the critical- wetting transition temperatures, 
but also the notoriously more difficult to predict first-order transition temperatures 
in the sequential-wetting scenario of hexane on brine. 

A crucial step in the description of the first-order transition temperature is the 
determination of the thickness 6 of the depletion layer, a thin film of pure water that 
forms near the brine/alkane interface. The two attempts made up to now to describe 
the wetting temperatures in the sequential wetting of hexane on brine rely on a multi- 
layer (four or five layers, respectively) DLP-type calculation involving layers of brine, 
water, dense liquid hexane in our case, liquid hexane, and vapor. In the first approach 
by Bertrand et al.^^ the free energy per unit area of such a four-layer structure 
is calculated, identified with the shift in contact energy caused by adding the salt, 
and subsequently converted into a shift of the first-order transition temperature as 
compared to the first-order wetting transition temperature of hexane on pure water 
in a conventional Cahn theory using a phase-portrait technique. The values of Tjjj^i 
resulting from this approach were shown to depend sensitively on the value of 
Excellent agreement with the experimental data was found for 5 ~ 2 A. 

Similarly, the predictions of T^^i within our approach presented in Sec. II. are 
in very good agreement with the experimental results if we choose 6 to be approx- 
imately 1.9 A. Interestingly, our theory is as sensitive towards changes in 6 as the 
one by Bertrand et al. (cf. Fig. 4). Therefore, both treatments have descriptive, 
but hardly any predictive power. We can conclude that the thickness of the layer of 
pure water is 1.8 - 2 A, but we would not be able to predict this thickness to within 
comparable accuracy for a different salt. In particular, there is a priori no reason 
why, for a different salt, the graphs of T^„,i(cNaCi) and T^,c(cNaCi) should be paral- 
lel; the two lines might actually meet in a critical endpoint. In view of the missing 
predictive power of the approach, i.e. without knowing 5 beforehand, however, there 
is no point in calculating the two lines of transition temperatures for another salt 
at the current state of theory. It is worthwhile noting that the value of 1.8 - 2 A 
for the hydrodynamic radius of Na^ and Cl~ is consistent with values determined 
in alternative way^ and reasonably well suited for describing the surface tension of 
NaCl solutions within the theory of Levin and Flores-Mena (cf. Fig. 3). 
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The major differences between the two approaches presented up to now, by 
Bertrand et al. and by us, are, firstly, in the way in which the changes due to 
the addition of sah are attributed to short-range and long-range forces, respectively, 
and, more importantly, in the reasoning why the depletion-layer thickness 5 should 
be constant. With respect to the former, our approach is to leave the contact energy 
(clearly a short-range interaction) unaltered as compared to the case of hexane on 
pure water since, as before, due to the presence of the depletion layer of thickness 5, 
only pure water is in direct contact with the adsorbate. We also retain the long-range 
forces between the layer of pure water and hexane; however, there are new terms in 
the long-range field which describe the interaction between the brine phase and the 
adsorbate across the 'layer' of pure water of thickness 5. All contributions to the 
free-energy functional in Eq. (jH^ are calculated at the actually relevant tempera- 
ture. Thus, our approach enables us to obtain absolute values of c and i if 5 is 
determined from independent measurements of the hydrodynamic radius or deduced 
by matching to the experimental (wetting) transition temperatures. 

The second and main difference of our treatment to the one proposed by Bertrand 
et al. is the justification of having a constant, i.e. salt-concentration independent, 
depletion-layer thickness 5. Bertrand et al. attribute the existence of this layer 
to electrostatic image charges that repel the ions from the brine/alkane interface. 
Based on this mechanism, Onsager and Samara^ had developed their theory of the 
increase of the surface tension of water on the addition of salt. They succeeded in 
deriving the correct limiting law for low salt concentrations and - as a by-product of 
their theory - gave an expression that relates the salt concentration to an equivalent 
depletion-layer thickness 5. The original Onsager-Samaras theory actually calculates 
the concentration profile and predicts an exponentially rapid approach of the local salt 
concentration to its bulk value as one moves away from the interface. By integration 
over this profile, one can compute the deficiency of ions near the interface. This value 
can, then, on the assumption that there is a layer which is completely devoid of ions 
(but, according to Onsager and Samaras, this should not be the case), be converted 
into an equivalent thickness of such a layer. It is this layer thickness that Bertrand 
et al. identify with 5 in their approach. Within Onsager-Samaras theory, however, 
this thickness is predicted to decrease slowly (logarithmically) but significantly with 
increasing salt concentration, so it is clear that the layer thickness cannot remain 
constant in the approach of Bertrand et al^ It is noteworthy that the equivalent 
depletion-layer thickness is only a by-product of the Onsager-Samaras theory and 
that it is not actually used to compute the excess surface tension of the electrolyte.!^ 

Our approach, in contrast, is based on the recent theory of Levin and Flores- 
Mena who argue that at higher salt concentrations (cNaCi > 0.15 mol/L and, there- 
fore, clearly in the range of salinities present in the sequential-wetting experiments) 
Onsager-Samaras theory neglects the 'intrinsic' depletion layer that is formed by the 
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hydration sphere of the ions (at low concentrations, this contribution is neghgible 
because the electrostatic repulsion creates a much thicker layer). In this picture, it 
becomes much more transparent why 6 is independent of the salt concentration (as 
long as there are enough water molecules for each ion to have a complete hydration 
sphere) and about 2 A in size. The local salt concentration will increase very rapidly 
to its bulk value in the region beyond the hydration-determined depletion layer thick- 
ness because electrostatic shielding is extremely efficient at high salt concentrations. 

In conclusion, the theory presented in this article - as well as the one by Bertrand 
et al^ - is able to describe the first-order transition temperature in the sequential 
wetting of hexane on brine if the thickness of the depletion layer is chosen to be about 
1.9 A. The predictive capacity of both approaches is, unfortunately, much lower. If 
one had a more precise a priori knowledge of the hydrodynamic radii of the ions, 
Tyj^i and T^, c could be predicted for solutions of different salts, which would be a 
major advantage because the actual experiments are very time-consuming.^ Such 
predictions might help identify suitable candidates of solutes for producing a critical 
endpoint in the wetting phase diagram, which is of particular interest .^^'^^'^ 
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TABLE I: Representative equations for the dielectric properties of water, brine, liquid 
hexane, and vapor. In addition to the relative static dielectric permittivity £:(0) and 
the refractive index ra, a characteristic absorption frequency common to all four media 
of Vf. = C(Je/(27r) = 3 X 10^^ s~^, which is in the UV range, has been used. T denotes 
the absolute temperature and CNaCi the salt concentration in the brine phase. 

medium relative static dielectric permittivity e{0) 

water=?9- 249.21 - 0.79069 K"! T + 0.72997 x lO^^ R-^ 

brin(P £:watcr(0) - 11 L mol"^ CNaCi 

hexaneP 1.890 - 0.00155 K'^ (T - 293.15 K) 

vapor ^ 1 

medium refractive index n 

wateiED 1.33436 - 1.50585 x 10"^ K^^ (T - 273.15 K) - 1.94586 x 10"^ K^^ 

X (T - 273.15 K)2 + 5.23889 x 10"^ K'^ (T - 273.15 K)^ 

brin^ riwator + 0.00918 L mol"^ CNaCi 

hexan(P ^z^^ter + 0.049 - 0.0004 K'^ (T - 273.15 K) 

vapor ^ 1 



TABLE II: Representative equations for the coefficients and 04, which are given 
by ttg = —W'/{6npi) and a'^ = B'a/^Anpi) as calculated from Eqs. (jH^ and (j^ . 
respectively, using the data listed in Table I, for various concentrations of salt. W and 
B' are the quantities actually listed in the table; W is represented in the form W = 
aw' + bw' {T — 273.15 K) + cw' (T — 273.15 K)^, where T is the absolute temperature, 
and B' is simply given by i?' = 2 W. 

CNaCi [mol L-1] W [J] 





aw [J] 




bw [J K- 


'] 


Cw [J K 






















0.5 


-5.064 X 10" 


-22 


5.329 X 10 


-25 


1.559 X 10 


-27 


1.0 


-1.012 X 10- 


-21 


1.069 X 10 


-24 


3.207 X 10 


-27 


1.5 


-1.516 X 10" 


-21 


1.609 X 10" 


-24 


4.998 X 10 


-27 


2.0 


-2.018 X 10- 


-21 


2.150 X 10" 


-24 


7.033 X 10 


-27 


2.5 


-2.519 X 10- 


-21 


2.689 X 10" 


-24 


9.509 X 10 


-27 
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FIG. 1: Part (a) shows the geometry of the four-layer structure water (1) /dense hquid 
hexane (4) /hquid hexane(3)/vapor(2) that is considered for the wetting of hexane on 
pure water. Part (b) extends this picture for the wetting of hexane on brine and 
introduces a fifth layer in a structure of the following type: brine(l) /water (4) /dense 
liquid hexane(5) /liquid hexane(3)/vapor(2). The film thicknesses of the various layers 
(4, 5, 3) are labeled 6, a, and I, respectively. 

100 
80 

^ 60 
O 

o 

^ 40 
20 


1 2 

FIG. 2: First-order and critical- wetting transition temperatures for hexane on brine as 
a function of the salt concentration. The open circles denote the experimental critical- 
wetting temperatures,^^ while the continuous line marks the theoretical prediction 
for this transition, as obtained from the temperature at which the Hamaker constant 
changes sign. The open diamonds represent the experimental first-order transition 
temperatures,'^ and the filled squares, connected by the dashed line to guide the eye, 
are our theoretical predictions using a depletion- layer thickness of (5 = 1.9 A. 
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FIG. 3: Excess surface tension of a 1-1 electrolyte, e.g. NaCl, over that of pure water 
as a function of the salt concentration. The filled circles represent the experimental 
data for NaCl at 25°C obtained by Matubayasi et al}'^ The lines are results of 
calculations based on the theory of Levin and Flores-MensP^ using different values 
for the hydrodynamic radius 5 of the ions; from top to bottom: 2.5 A (dash-dotted), 
2.125 A (thick long-dashed), 1.9 A (continuous), 1.75 A (thin dashed). 
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FIG. 4: Predictions of the first-order transition temperatures in the sequential wetting 
of hexane on brine for various values of the depletion-layer thickness 6; from top to 
bottom: 2.5 A (open circles, dash-dotted line), 2.125 A (filled triangles, thick long- 
dashed line), 1.9 A (filled squares, thin continuous line), 1.75 A (filled circles, thin 
dashed line), concentration-dependent 6 according to Onsager-Samaras theorjl^^*^ 
(open diamonds, thin long-dashed line). The thick continuous line marks the critical 
wetting transition temperatures as calculated from the Hamaker constant. 
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